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We examine an inflationary model in R + R^ gravity with torsion, where R^ denotes five 
independent quadratic curvature invariants; it turns out that only two free parameters remain 
in this model. We show that the behavior of the scale factor a{t) is determined by two scalar 
Q^, fields, axial torsion x(t) and the totally anti-symmetric curvature E{t), which satisfy two first- 

5^ ' order differential equations. Considering % « during inflation leads to a power-law inflation: 

5^1 a ~ (i + Ay where 1 < p < 2, and the constant A is determined by the initial values of E, 

X and the two parameters. After the end of inflation, x and E will enter into an oscillatory 
CO ■ phase. 

> 

\o 
o 
p 

o 

00 ' Submitted to: Class. Quantum Grav. 

O' 

> 

X 



PACS numbers: 04.50.-h, 04.50.Kd, 98.80.Cq, 98.80.Jk 



Inflation in R + R"^ Gravity with Torsion 2 

1. Introduction 

Since cosmic microwave background radiation (CMB) was discovered, the standard hot big 
bang cosmological model becomes widely accepted. However, the hot big bang Universe 
leads to several problems, e.g. horizon and flatness problems [0|. Inflationary scenario was 
then proposed to resolve these problems (see ifTOl '). The standard approach to inflationary 
Universe is to introduce scalar inflaton fields without modifying general relativity (GR). 
However, these inflaton fields will violate the strong energy condition. Instead of introducing 
inflatons, it has been found that several alternative theories of gravity can also provide 
inflationary Universe. A well-known pure gravity inflationary model is based on a fourth- 
order gravity theory TZ + eR?, i.e., adding a quadratic pseudo-Riemannian scalar curvature 
'R? to the Einstein-Hibert action without introducing any scalar inflaton field [f2Tl l23l [TSl . 
In the following, we put hat on any geometric object W to indicate that W is defined from 
Levi-Civita connection V. In [|23l (see [[22l| for further generalization), it provides a better 
understanding of IZ^ inflation by showing that the Lagrangian — 2A + 1Z + eR? is conformally 
related to Einstein gravity plus a scalar field, where A denotes the cosmological constant. 
Furthermore, it has been demonstrated that the Lagrangian — 2A + 7?. + eR? can be considered 
as a general quadratic curvature Lagrangian for homogeneous, isotropic cosmological model, 
i.e., pseudo-Riemannian curvature square iZ°'^^'^iZabcd and Ricci curvature square iZ°'^iZab can 
be expressed in terms of 7^^ by using the Gauss-Bonnet theorem and a homogeneous, isotropic 
assumption [2J. It should be stressed that this demonstration will not be true in Riemann- 
Cartan spacetime, i.e., metric-compatible connection Vg = with torsion. 

Since magnitudes of curvature are expected to be very large in the early Universe, the 
effects of quadratic curvature invariants may dominate. So it is reasonable to include them in 
a Lagrangian, which is a linear function of scalar curvautre. However, it is not satisfactory 
to restrict these investigations in pseudo-Riemannian geometry, since the connection V and 
the metric g are fundamentally independent. Therefore, the fundamental variables of gravity 
should be considered as {g,V}. It motivates us to investigate quadratic curvature effects, 
which is based on Riemann-Cartan spacetime, during inflationary epoch. In this paper, we 
extends the Einstein-Cartan theory to include all of quadratic curvature invariants, i.e. five 
independent invariants, in the early Universe. Since dimension of torsion and curvature are 
L~^ and L^^ respectively, quadratic torsion effects can be neglected by comparing them to 
quadratic curvature in the early Universe. So we do not include quadratic torsion terms in 
the Lagrangian. It turns out to be a degenerate case of the Poincare gauge theory of Gravity 
(PGT), which does have three quadratic torsion invariants [|6||. A homogeneous, isotropic 
cosmological model in the PGT has been developed and many special cases and solutions 
are discovered [81 [161 [l2l l20ll- These special cases have been used to explain dark energy 
problems [[T7ll20l l and inflation fT6\. In fT6l, the authors reduced the PGT field equations 
into a second-order differential equations with respect to Hubble radius by certain restrictions 
on indefinite parameters and also introduced a scalar field and matter fields during inflation. 
However, our inflationary model is completely different from ref [fT6ll . We start from a general 
quadratic curvature Lagrangian in the Riemann-Cartan spacetime and then obtain a pure 
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gravity inflationary model without introducing any scalar or matter fields. 

Our purpose in this paper is to investigate quadratic curvature effects in the early Universe 
and also to analyze whether these effects will generate inflation. Since our analysis is based 
on the Riemann-Cartan spacetime, there are five independent quadratic curvature invariants 
in the Lagrangian. It is worth to point out that our work is not a special case of f(JZ) gravity 
with torsion [4], where 7^ is a scalar curvature. In the f{TZ) gravity with torsion, it does not 
include all possible curvature invariants, e.g. Ricci square TZ^^TZab, in its Lagrangian. There 
is no fundamental reason for us to neglect these terms in the early Universe. Moreover, we 
do find that these neglected terms can yield a power-law inflation. The field equations of 
f{TZ) gravity with torsion give totally antisymmetric torsion tensor, i.e. axial-vector torsion, 
vanishing. Conversely, our inflationary solution is purely determined by the time-component 
of axial- vector torsion x called axial torsion. 

The plan of this paper is as follows. In Section 2, we give a brief review of Riemann- 
Cartan geometry. Moreover, we start from a general quadratic curvature Lagrangian and 
obtain an effective quadratic curvature Lagrangian, which contain only three quadratic 
curvature invariants, for a homogeneous, isotropic cosmological model. A similar argument 
has been done in pseudo-Riemannian spacetime [|2l[l5l. In Section 3, we derive the equations 
of motion for inflationary cosmology. Since we assume that there is no matter field during 
inflation, it further reduces one more parameter. The equations of motion will contain only 
two parameters. In Section 4, we find an inflationary solution which are corresponding to a 
power-law inflation. The e-folding number N, which is determined by one of two parameters, 
can satisfy the standard requirement A^ > 60. At the end of inflation, x will enter into an 
oscillatory phase. In this paper, h = c = 1 and the metric signature is ( — h ++). 

2. Riemann-Cartan geometry and the quadratic curvature Lagrangian 

Since Einstein developed GR in terms of a pseudo-Riemannian spacetime, it was argued that 
GR has a pre-geometry, especially when one discovered that intrinsic spins can naturally 
generate torsion tensor. A directly extended theory of GR to include torsion is called Einsten- 
Cartan theory, which is based on a Riemann-Cartan spacetime. If the spin densities of matter 
fields vanish, the Einstein-Cartan theory will return to GR. Since torsion effects generated by 
spin sources is too small, the current experiments still cannot distinguish their differences. 
However, these two theories may become significantly different in early Universe due to 
quadratic curvature effects. 

The geometry of a Riemann-Cartan spacetime is described by metric g and metric- 
compatible connection with torsion V. The curvature 2-forms i?"^ and torsion 2-forms T" 
can be expressed in terms of connection 1 -forms uj°'i, and an orthonormal co-frame e": 

R\ = du\ + u\Au\ (1) 

T" = de" + a;";, A e^ = De^ (2) 

where D is the co variant exterior derivative |l3l. In an orthonormal frame, Uab satisfy 
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antisymmetric condition 

^ab + t^fea = 0, (3) 

where the indices are lowered by rjab = diag(— 1, 1,1,1). In the following, Latin indices are 
raised and lowered by r]ab- 

It has been shown that the curvature 2-forms Rab can be irreducibly decomposed into 

12 3 4 

the following pieces [143 : Rab =Rab (Weyl)+ Rab (Paircom)+ Rab (Pseudoscalar)+ Rab 

5 6 

(SymRicci)+ Rab (AntiRicci)+ Rab (Scalar). Moreover, the irreducible components of 

12 3 

torsion 2-forms yield T"" =T'' (Tensor) + T" (Vector) + T" (Axitor). It should be pointed 

2 3 5 

out that one has Rab=Rab=Rab= and T" = in the pseudo-Riemannian spacetime. In 

a 

terms of Rab, the most general quadratic curvature action (neglected any possible derivative 
on curvature) in the Riemann-Cartan spacetime is given by 

S[e'',u;\] = ^ [ [A*l + Co Rab A *e^' + R""' Aj^Ca* Rab], (4) 

where k = 8nG and e""^c-d denotes e'* A • ■ ■ A e^ A Cc A • ■ ■ A e^. A is the cosmological 
constant and * denotes the Hodge map associated with g. cq is a dimensionless parameter 
and [cq] = L^.ln the PGT, the most general gauge invariant action should also includes three 
irreducible pieces of torsion square in its Lagrangian |[6l[8l. Although Rab have six irreducible 
components, there are only five independent quadratic curvature invariants due to generalized 
Gauss-Bonnet theorem, which gives [[1911 

SabcdR"^ A R"^ = an exact form. (5) 

Using ([5]), one of irreducible components of Rab can be replaced in terms of an exact form 
which does not contribute to field equations. So it turns out to have only six parameters 
(including cq) in the field equations, which are obtained from varying dH). In order to obtain 
the Newtonian theory in weak field approximation ||5l [HI, we should fix cq = 1, which 
actually gives Einstein-Cartan action with A if quadratic curvature invariants vanish. 

It is straightforward to obtain field equations by varying S with respect to {e'*} and {cuab} 
and these equations should contain five independent parameters |j However, one can further 
reduce two parameters in the field equations according to an assumption of homogeneous, 
isotropic Universe. It can be verified that the homogeneous and isotropic assumption yields 

1 5 1 5 

Rab=Rab= by using (IS) and uvn . so quadratic Rab and Rab do not contribute to field 
equations for homogeneous, isotropic cosmology. It is worth to point out that they may 
participate in the perturbed field equations |§| Since our current work is only concentrated 
on homogeneous and isotropic Universe, we can start from the following effective action 

SEie", uj\] = ^ [ [A* 1+ Rab A *e'^^ + h Rab A e^' A *{Rab A e^')] 

+ b^PaA *P" - be Rab A *e^' A *{Rab A *e^')], (6) 

I These field equations are equivalent to the field equations (3.3) and (3.4) in [8] by neglecting torsion square 

terms and rearranging parameters of quadratic curvature invariants in (3.3)-(3.4). 

§ The cosmological perturbation analysis for R + R^ gravity will relate to our future works. 
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where Pf, = iaR°'b denotes Ricci 1 -forms. So it only involves three indefinite parameters. In 
the following, we set A = since the effects of A in the early Universe are expected to be 
small. 

Varying 5"^; + S'matter with respect to the orthonormal co-frames {e"} and connection 
1-forms [uab] gives following field equations: 

Rab A *e''\ - h E{E * Ce + 2 i.Rab A e'^') + h{ tc{P' A *Pb) 
+ 2 icR"'' Aia* Pb) + h n{2Rab A *e''\ ~n*ec) = -2k t^, (7) 

r A *€'"'' c + 2 b3D{E e"^) - b^D * {P^ A e" - P'^ A e'') 
+ 2 beD{n * e"'') = -2k 5°^ (8) 

where IZ = iaP"" is the scalar curvature and E = \Rabcd^°'^^'^ is the totally antisymmetric 
curvature, ia denotes the interior derivative and e"^'^'^ is the Levi-Civita antisymmetric e- 
symbol. Tc and 5*"^ are stress-energy and spin density 3-forms of the matter fields. 

3. Equations of motion for inflationary cosmology 

In an isotropic, homogeneous cosmological model, the spacetime metric and torsion 2-forms 
become JH: 

3 

9 

A=l 

T" = /(t)e"Ae° + 2x(t)*(e'^Ae°), (10) 

where curvature constant fc= -1, 0, or 1. For simplicity, we only concentrate on the flat 
Universe A; = 0. Moreover, the non-vanishing components of r^ and 5""^ yield 

Ta ='nOaipit)+pit))*eo+pit)*ea, (11) 

S^^ = - q{t) e^^ A *(e''l A e°) + s{t) e° A e" A e^ (12) 

where square bracket indicates index anti-symmetrisation. In standard inflation models, the 
energy density p and pressure p of an inflaton field are given by p = ^0^ + F(0) and 
p = \(tP' — V {(f)) ^. It still lacks a fundamental theory which can explain the origin of 
the inflaton fields and also derive the potential V{(i)). Moreover, there is no observational 
evidence to yield any feature of matter fields during inflationary epoch. In this paper, we will 
assume r^ = 3°"^ = in the early Universe and then obtain a pure gravity inflationary model 
in the Riemann-Cartan spacetime. 

Using ©-([111), the field equations ©-dH]) with A; = yiel(i||] 

{{H + ff - x'} + |i?' - ^b,Ex{H + /) 

-^^n^^^-^^^^n{iH + ff-x^} = ^-^, (13) 

n = K{p-?>p), (14) 

II (fT3Tl-(fT6ll is a degenerated case of the field equations (4.2)-(4.5) in ||8]. 



1 ^ 

dt ® dt + a^{t) (1 + -kr'^Y^ XI ^^^ ® ^^^' ^^^ 
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{be + |)(^ - 27^/) + 2hx{Hx + x) + ^hEx - / = - y, (15) 



^4 + 666;^,, or, ,,r/rr , f\2 ,21 



hiE - 2fE) - ^ , °^ 7^x - 2b4X{(^ + fr -X'}~X = -^^s, (16) 



where 



n = 6iH + 2H^ + 3Hf + f + f- x'), (17) 

E = Qix + ^Hx + 2xf), (18) 

and H = - denotes the Hubble parameter. It can be shown that (fT3])-(fT6l) form a complete 
system to describe the evolution of the isotropic, homogeneous flat Universe once the equation 
of state p = p{p) is given. In particular, p = p = (inflationary epoch) or p = ^ (radiation 
domination) both yield 7?. = (see (fl?)) ). It is not difficult to see that (fT3l) yields the Friedmann 
equations when the quadratic curvature terms and the spin density vanish. By substituting 
7?. = into (fT3l) . (fT5l) and (fT6l) . one can easily verify that only ^3 and 64 left in these equations. 
Since we consider p = p = q = s = during inflation, (fT3l)- (fT6l) with (fTSi) then become 

E = ^ + 2{<!> - H)E + ^X{'^' - X'), (19) 

03 03 

X=^-Hx-2x^, (20) 

o 

with two algebraic equations: 

^=2bsEx±A, (21) 

H= - |ex + 863 64 £^x' ± (1 + 4fo4x')^, (22) 



where ^ = Ji^b^ExY + x^ - ^ and $ = iJ + /. We further simplify ([l9l) and ^ by 
substituting (|2TI) and (l22l) into them, which gives 



E = ^ + (463 + ^) E\, (23) 

63 3 

^ = ? + (t - 463)i?x' T (3 + 464X')^X - 863ft4i?x', (24) 

It turns out that the evolution of the inflationary Universe is characterized by 63, 64 and the 
initial values (t = 0) of -E and x- We should point out that here t = denotes the beginning of 
inflation. It is reasonable to require 63 = — 6 < 0|5llin order to ensure that ^ is a real function. 
Moreover, we assume that the initial values Eq and xo are both less than or on the order of the 
Planck scale to ensure the classical validity of the evolution [fTSl . 

4. Power-law inflation 

We first notice that (|22|)- (|24|) has a similarity to hybrid inflation models [fTTl [T2|, with the 
effective potential given hy V{E,x) = Sk^^H^. In hybrid inflation models, the effective 
potentials V{a, 4>) are constructed in such a way that one scalar field a = during slow-roll 
inflation, and starts to roll down to its true vacuum when the other scalar field (j) falls to a 

^ The requirement of positive kinetic energy in the spin 0^ mode of Hnearized PGT also gives 63 < ID. 
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critical value 0c. We adopt a similar arguement by considering |x| ~ |xo| ^ Ch during 
inflation and I-EqI ~ Ch ■ By neglecting x^ in A, (l24l) reduces to an algebraic cubic equation 

(32„ + -)y'-(20+---)x'^-(- + -)X + -«0. (25) 

where x = ^Xo > ^^^ '^ — ~^' which are both dimensionless parameters. From 
numerically and analytically analyzing (l25l) . we find that the condition a ~ x ~ C)(l), which 
is consistent with neglecting x^ in A, can naturally provide the sufficient number of e-folding 
(see below), so we will assume the condition being satisfying in the following discussions. 
Since 6 = ^ ^ /pj^, one can neglect f- in (|23l) . and it becomes 

E^-{4+^)bxoE'. (26) 

(|26l ) then yields a general solution 

B«b(* + |-)". (27) 

where i3 = (^2+1) b ■ ^^ substituting (|27l) into (l22l) . we obtain // :^ p (t + J-)"^, where 



for the case xoE > 0, and ± should be changed to =F for the other case xoE < 0. So a 
solution of a(t) for a(t) > gives 

a^Aoit + ^y, (29) 

Eq 

where Aq = ao(;f-)^^ > 0, and oq > denotes the initial value of a. 

It is known that a general criterion for inflation is d > 0, which may be sufficient to solve 
the horizon and flatness problems. It turns out that to require p > 1 in (|29|) corresponds a new 
type of the power-law inflation [HI [HI. From (l28l) . one can see that p > I will limit the value 
of a. For example, if we consider a degenerate case a = 0, it gives p = | which does not have 
inflation. However, in the a = 1 case, which yields x ~ 0.623, the numerical calculations 
of (l28l) while choosing the negative sign in the round bracket yields p ~ 1.938, which gives 
power-law inflation. From numerical calculations, we also find that I < p < 2 when a > 0. 
In addition to d > 0, the amount of inflation is usually required to satisfy A^ > 60, where 
]\f = ln ^iissdi denotes the number of e-foldings, and tend means the end of inflation. In order 
to estimate A^, we consider that the inflation comes to an end when (|25]) is no longer valid, 
and it occurs when x^ in A cannot be neglected , i.e. 

IbElJr^xl or \b'El^\r^Oil). (30) 

From (l27l) and (l30l) with the help of a; ~ x ~ 0(1), the values i^cnd and tend can be estimated 
as Ecnd ~ (i2+a)b ^^^ ^cnd ~ — . By Substituting tend ~ — into the number of e-foldings 

iV = ln('^tend + l] , (31) 
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we obtain e^ ~ ( (4+ |)£'o b)^- It turns out that the value of A^ is determined by the parameter 
b. In the a = 1 case, A^ > 60 requires b > 10^^ l^^ ~ 10~^^ m^, which satisfies the tests of 
solar experiments [7J. After the end of inflation, x and E both become very small, and the 
linear terms in (fT9l) and (l20l) should finally dominate. In the linearized equations of (|23l) and 
(|24)) . one can easily obtain oscillatory solutions: x ~ sin -4t and £' ~ cos -^. Whether 
these oscillatory phases can provide a reheating process is still in progress. 

5. Conclusion 

In conclusion, we first discovered a pure gravity inflationary model in i? + i?^ gravity theories 
with torsion, and presented a detail analytic analysis of this model. It turns out that the 
dynamics of the inflationary model is purely determined by two parameters b and a with 
two initial values Eq and xo- We further found that |xo| ^ Ch and \Eq\ ~ /pj^ yields a 
power-law inflation if « > 0. In this power-law inflation solution, x is near a constant xo and 
E ~ t^^. Since x cannot be near xo after E decays to satisfy the condition \b E^^^^\ ~ Xq, 
it can be considered as the end of inflation. The period of inflation can be estimated to be 
— . Moreover, the e-folding number A^ is completely determined by b and also satisfy the 
requirement A^ > 60 when b > 10^^ l^^i ~ 10~^^ m^. 

After the end of inflation, the system may enter into an asymptotical regiem which yields 
an oscillatory phases. These oscillatory phases are expected to give a reheating process. Since 
CMB anisotropy has been discovered, a further understanding of the reheating process and 
primordial perturbations will provide more restricted constraints on b, a, Eq and xo- These 
issues will be considered as our future work. 
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